Phase separation in coupled chaotic maps on fractal networks 
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The phase ordering dynamics of coupled chaotic maps on fractal networks are investigated. The 
statistical properties of the systems are characterized by means of the persistence probability of 
equivalent spin variables that define the phases. The persistence saturates and phase domains freeze 
for all values of the coupling parameter as a consequence of the fractal structure of the networks, in 
contrast to the phase transition behavior previously observed in regular Euclidean lattices. Several 
discontinuities and other features found in the saturation persistence curve as a function of the 
coupling are explained in terms of changes of stability of local phase configurations on the fractals. 

PACS numbers: 05.45.+b, 89.75.Kd 
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Coupled map lattices have provided fruitful and com- 
putationally efficient models for the study of a variety 
of dynamical processes in spatially distributed systems 
. The discrete-space character of coupled map systems 
makes them specially appropriate for the investigation 
of spatiotemporal dynamics on nonuniform or complex 
networks. Phenomena such as pattern formation, spa- 
tiotemporal intermittency, nontrivial collective behavior, 
synchronization, etc., have been extensively studied in 
coupled map systems defined on fractal lattices hier- 
archical structures trees Q, random graphs |j|7 small- 
world networks HJ, and scale- free networks Q. 

Recently, there has been much interest in the study 
of the phase-ordering properties of systems of coupled 
chaotic maps and their relationship with Ising models in 
statistical physics 0, 0, 0, 0, 0] . These works have 
invariably assumed the phase competition dynamics tak- 
ing place on a uniform Euclidean space; however, in many 
physical situations the medium that supports the dynam- 
ics can be nonuniform on some length scales. The nonuni- 
formity may be due to the intrinsic heterogeneous nature 
of the substratum such as porous or fractured media, or it 
may arise from random fluctuations in the medium. This 
article investigates the phenomenon of phase ordering in 
coupled chaotic maps on fractal networks as a model for 
studying this phenomenon on nonuniform media. The 
class of fractal networks being considered corresponds to 
generalized Sierpinski gaskets (GSG) embedded in Eu- 
clidean spaces of arbitrary dimension d 2\. In particular, 
this model of coupled maps on fractal networks yields a 
situation to explore the role that the connectivity of the 
underlying lattice plays on the statistical properties of 
phase ordering processes in nonlinear coupled systems. 

Deterministic fractal networks, such as GSG, can be 
generated in any d-dimensional Euclidean space as fol- 
lows 0- At the nth level of construction, the fractal 
consists of N = (d + l) n d-dimcnsional hypertetrahe- 
dral cells whose coordinates can be specified by a se- 
quence (a\a.2 ■ ■ ■ a n ), where a m can take any value in a 
set of (d+1) different symbols which can be chosen to be 



{0,1, ... ,d}. At level n+1, each cell (aia.2 ■ ■ ■ a n ) splits 
into d+1 cells scaled down by a longitudinal factor of two, 
and which are now labeled by (a\a2 ■ ■ ■ a„a n +i), where 
the first n symbols of the sequence are the same as the 
parent cell. Given this construction rule, the fractal di- 
mension of the GSG is dt = \og(d + l)/log2. A label 
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(cvia2 ■ ■ - a n ) can be written as (a±a2 ■ ■ ■ a„ 
for some s € {1,2, ...,n}, where af means the se- 
quence of s identical symbols a.;. The cell with this la- 



bel has a neighborhood Afr ai 



(aiCK2 ■ 



with d+1 elements 



labeled by (evi . . . a n -i(a n + j)) (j — 1,2, ...,d) and 
(a i . . . a n - s +ia s n _ s ), where the addition ai+j is defined 
modulo (d+1). If s — n, then the cell is one of the d+1 
vertices of the gasket, labeled by (a™), and it has only d 
neighbors belonging to the same parent cell. An integer 
index i = 0, 1, . . . , (d + 1)" — 1, can be assigned to each 
cell of the lattice at the level of construction n by the 
™\ei = Yl=iUm{d+l) n - m . 

The equations describing the dynamics of the dif- 
fusively coupled map system defined on these fractal 
networks at level of construction n, embedded in a 
d-dimensional Euclidean space, are 



x f+ i(ai . . . a„) = (1 - e)/ (x t (ai . . . a n )) + 

3TT S(/3 1 ... / 3„)G7V (Q1 ... Q , l) / ( X t(Pl ■ --AO) i 
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where Xt{ot\ . . . a n ) gives the state of the cell (ai . . . a n ) 
at discrete time t; (pt\ . . . a n ) and (/3j . . . j3 n ) label the 
(d+1)™ cells on the gasket; e is a parameter measuring the 
coupling strength between neighboring sites, and f(x) is 
a nonlinear function that expresses the local dynamics. 
Equation (JJJ also applies to the (d+1) vertex cells of the 
fractal network, except that the coefficient of the sum is 
e/d, since each of these cells have d neighbors. 

As local dynamics, we assume a piecewise linear, odd 
map Q 



f(x) = 




1,-1/3] 
1/3,1/3] 
[1/3,1]. 



(2) 



2 



When the parameter n <E [1,2], the map possesses two 
symmetric chaotic attractors contained in the intervals 
I ± = [±/x(2 — /x)/3, ±/x/3], and separated by a gap 
7° = [— /i(2 — fi)/3,n(2 — A*)/3]. For values of /z close 
to two, the size of the chaotic intervals is larger than the 
gap. Then the local states have two well defined sym- 
metric phases that can be characterized by spin variables 
defined as the sign of the state at time t, at{a\ . . . a n ) = 
sign(a; t (ai . . . a„)). 

To study the phase-ordering phenomenon of the cou- 
pled maps on fractal networks we fix the local map pa- 
rameter at the value = 1.9 and set the initial condition 
as follows: if the number of cells (d + 1)™ in a lattice is 
even (d odd), exactly one half of the sites are randomly 
chosen and assigned random values uniformly distributed 
on the interval I + , while the other half are similarly as- 
signed values on I~ . If the number of cells in a lattice 
is odd [d even), then the state of the remaining cell is 
assigned at random on cither interval I + or I~ . 

The statistical properties of the phase-ordering process 
on the fractal networks can be characterized by using the 
persistence probability Pt, defined as the fraction of cells 
which have not changed sign up to time t |l3j |. Figure 1 
shows Pt as a function of time for the GSG embedded in 
d = 3, for several values of the coupling parameter. For 
some ranges of the coupling, the persistence saturates 
in a few iterations, while for some other ranges of e, P t 
reaches its saturation value more slowly. In contrast, in 
regular Euclidean lattices the persistence saturates for 
small couplings, while it decays algebraically in time for 
coupling strengths greater than some critical value Q . 
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FIG. 1: Persistence probability as a function of time for the 
Sierpinski gasket embedded in Euclidean dimension d = 3 for 
different values of the coupling e. Lattice size is TV = 4 9 . 
Dotted curves correspond to values of e for which P t reaches 
its saturation value quickly. 

The domains formed by the two phases on fractal lat- 
tices reach a frozen configuration for all values of the cou- 
pling e £ [0,1]. Figure 2 shows the asymptotic patterns 
of the phase separation process on a GSG embedded in 



d = 3 at level of construction n = 3, for several val- 
ues of the coupling. Note that the configuration of the 
blocked phase domains changes as the coupling is varied. 
The domain configurations can be characterized by the 
fraction of sites in a given phase that have k neighbors 
in that same phase at time t, denoted by F t (k), with 
k = 0, 1, . . . , d + 1. For example, consider the pattern 
displayed in Fig. 2(a) where there are several sites in one 
phase, as those indicated by arrows, having all of their 
four neighbors in the opposite phase. Thus the asymp- 
totic fraction Poo(0) is greater than zero in this case. In 
Fig. 2(b), as e increases, Poo(0) becomes zero, but Poo(l) 
is finite since there are sites in one phase, as those sig- 
naled by arrows, having just one neighbor in that same 
phase. 




FIG. 2: Phase separation on the GSG embedded in Euclidean 
dimension d = 3 at level of construction n = 3, for several 
values of the coupling. Dark (light) color corresponds to the 
positive (negative) phase. Arrows signal local configurations 
described in the text, (a) e = 0.20; (b) e = 0.26; (c) e = 0.36; 
(d) e = 0.72. 

The relationship between the asymptotic behavior of 
the persistence and the frozen domain configurations on 
fractals becomes manifest in Fig. 3(a), which shows the 
saturation value of the persistence, P^, as well as the 
fractions Poo(0), Poo(l) and F 00 (2), as functions of the 
coupling parameter for the GSG embedded in Euclidean 
dimension d = 3. Several discontinuities are observed 
in the curve of Poo in Fig. 3(a). The first discontinu- 
ity of Pqo occurs at the value of the coupling e Q = 0.24 
where the fraction Foo(0) vanishes. This implies that lo- 
cal blocked configurations where a site has no neighbors 
in its same phase, as those indicated in Fig. 2(a), be- 
come unstable at the value of coupling e Q . The second 
discontinuity of Poo takes place at the value of coupling 
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continuities on each curve are related to the loss of sta- 
bility of the configurations where one local element has 
a majority of its neighbors in the opposite phase. The 
number of those configurations is J = nint((d + l)/2), 
where nint(a;) rounds x to its nearest integer. These J 
discontinuities take place at increasing values of the cou- 
pling ek for which the asymptotic fractions FooQt) van- 
ish, with k < J. The first discontinuity is related to the 
vanishing of the fraction Foo(0) at the value e — 0.24, 
independently of the embedding dimension d. This in- 
dependence is due to the presence of the normalization 
factor (d + in the coupling term in Eq. Q for all 
the lattices. For embedding dimensions d odd, there oc- 
cur a minimum of when the fraction Foo(J) becomes 
zero. In those cases the local configurations losing stabil- 
ity are those consisting of a site with half of its neighbors 
in one phase and the other half in the opposite phase. 
These local configurations are symmetric, in the sense 
that a change of phase of that site does not alter the 
phase composition of its neighborhood. These symmet- 
ric configurations can not happen in GSG associated to 
d even, and for those networks the minimum of Poo co- 
incides with the the last discontinuity that takes place 
when F^J — 1) vanishes. 

Note that for each fractal lattice there is a value of the 
coupling e max where a maximum of Poo is observed. The 
origin of such maximum can be analyzed through the 
average fraction of neighboring pairs that have opposite 
phases at time i, defined as 
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FIG. 3: (a) Left scale: Saturation persistence as a func- 
tion of e for the GSG embedded in Euclidean dimension d — 3, 
represented by the thick dark line. The values e , 61,62, and 
E mol are indicated by dotted vertical lines. Right scale: Frac- 
tion F t (k) vs. e. Circles Poo(O); triangles Fao(l); squares 
F aa (2). Lattice size is N = 4 9 . (b) as a function of e for 
GSG embedded in different Euclidean dimensions d. Contin- 
uous line: d = 2, N = 3 1 ; dotted line line: d — 5, TV = 6 7 ; 
dashed line: d = 7, N = 8 6 ; circle line: d = 11, N = 12 5 . 
The first discontinuity at e = 0.24, common to all lattices, is 
indicated. 



€\ = 0.30 where -Foo(l) becomes zero, and it is related to 
the loss of stability of local frozen domain configurations 
as those signaled in Fig. 2(b). In addition, reaches 
a minimum at the value of coupling £2 = 0.36, where the 
fraction Foo(2) decays to zero. For e > 62 the domains 
of the two phases grow in size as seen in Fig. 2(c) and 
2(d). The phase domains also form faster, reducing the 
number of phase switching of the elements and therefore 
producing an increment in the saturation values Poo up 
to a maximum occurring at e max = 0.72. 

A similar behavior is observed for all the fractal net- 
works embedded in different Euclidean dimensions d. 
Figure 3(b) shows the saturation persistence Poo as a 
function of e for GSG associated to different d. The dis- 
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where S(a t (i),a t (j)) = 1 if a t (i) = a t {j), and 
5((Jt{i),(Jt(j)) = otherwise. In Fig. 4 we show G\ 
as a function of the parameter e for different fractal lat- 
tices. The maximum of G\ for each lattice takes place at 
the value e max at which the corresponding curve of Poo 
reaches a maximum. When G\ is maximum the proba- 
bility that there exist cells in one phase having all of their 
neighbors in the opposite phase at the first iteration is 
also maximum. Therefore, for the value of coupling e max 
there is a greater chance that at the next iteration such 
cells change phase instead of their neighbors, yielding 
more stable domains with fewer elements of the network 
having to switch their initial phase. Consequently, the 
persistence probability, that measures the number of el- 
ements that have not changed phases, is maximum at 
that time. Since the domains that are being formed are 
the most stable, in successive times the persistence will 
sustain a maximum value for the value of coupling e max 
corresponding to each lattice. For e > e max the coupling 
is strong enough to induce transient changes in the phase 
of elements having the majority of their neighbors in that 
same phase and therefore producing lower values of Poo- 
The local effect captured by the quantity Gi in frac- 
tal lattices also appears in regular Euclidean lattices, 
although the asymptotic behavior of the persistence is 
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FIG. 4: Gi as a function of e for GSG embedded in different 
d. Dashed line: d = 2; continuous line: d = 3; dotted line: 
d — 5; circle line: d = 11. Sizes are the same as in Fig. 3. 
Dotted-dashed line represents Gi for a two-dimensional reg- 
ular lattice, (2D). 



different in those two network topologies. Figure 4 in- 
cludes the calculation of G\ as a function of e for a two- 
dimensional regular lattice; the maximum of G\ in this 
case occurs at e max = 0.67. This is the critical value 
of the coupling parameter found in Ref. Q, after proper 
normalization, for the phase ordering transition in the 
scaling behavior of the persistence in a two dimensional 
Euclidean lattice. At e = 0.67 the blocked states in the 
two-dimensional regular lattice give place to growing do- 
mains of the two phases, separated by a continuous inter- 
face. The interface motion is driven by curvature effects 
that cause changes of phase in many elements of the sys- 
tem and therefore a temporal decay in the persistence 
probability jij. In general, in regular Euclidean lattices 
the ratio between the length of the interface to the size 
of a domain decays as r^ 1 , where r is the average radius 
of the domain. On the other hand, in the fractal net- 
works the interface consists of a few disconnected cells 



separating large domains and no curvature can be de- 
fined. Because of the self-similarity of the structure the 
number of elements in a domain grows as r d f = [d + 1) , 
where / is an integer, while the size of the interface is of 
the order of (d + 1), independently of the domain sizes. 
Thus the ratio of the interface to domain size decreases 
as r -df(l-i)/l m foe fractals. This decay is much faster 
than in regular Euclidean lattices and accelerates with 
increasing domain size, forming stable separated phase 
domains. As a consequence, domains always freeze on 
the fractal networks and the phase transition observed 
in the temporal behavior of the persistence in regular 
Euclidean lattices does not occur in fractals. 

In summary, we have found that phase domains in 
chaotic maps coupled on fractal networks always reach 
a frozen configuration, causing the saturation of the per- 
sistence in time for all values of the coupling parameter, 
in contrast to the phase-transition behavior of the persis- 
tence observed in Euclidean regular lattices. The fractal 
nature of the spatial support is also reflected in the dis- 
continuities observed in the Poo vs. e curves in Fig. 3. 
The phase configurations of the local neighborhoods have 
similar transient manifestations in fractal networks and 
in Euclidean regular lattices, as seen in the emergence of 
a maximum of G\ at a value of coupling e max . However, 
the asymptotic and global properties of the phase order- 
ing process on these two network topologies are quite 
different, even when the number of local connections in 
the neighborhood is the same, as it happens for the two- 
dimensional Euclidean lattice and the GSG embedded 
in d = 3. The necessary correlations between elements 
for building complex collective dynamics are more likely 
to occur in the denser Euclidean lattices. These results 
suggest that the topology of the network and not the 
number of local connections or the dimensionality of the 
space determine the asymptotic collective behaviors that 
may emerge on networks of coupled chaotic elements. 
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